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A theoretical study of tiie subpicosecond fragmentation of Ceo clusters in response to ultrafast 
laser pulses is presented. We simulate the laser excitation and the consequent nonequilibrium relax- 
ation dynamics of the electronic and nuclear degrees of freedom. The first stages of the nonequilib- 
rium dynamics are dominated by a coherent breathing mode followed by the cold ejection of single 
C atoms, in contrast to the dimer emission which characterizes the thermal relaxation. We also 
determine the nonequilibrium damage thresholds as a function of the pulse duration. 



During the last decade much attention has been paid to 
the study of the ultrafast fragmentation of molecules and 
clusters upon excitation with femtosecond laser pulses 
The sub-picosecond relaxation dynamics of the 
nuclear degrees of freedom, as a typical ultrafast phe- 
nomenon, becomes interesting when it presents new fea- 
tures, qualitatively different from the well known thermal 
relaxation processes. There are many recent examples 
of such new phenomena in solids, like the cold melting 
of graphite and silicon |^,^, ultrafast desorption ||^ or 
femtosecond graphitization of diamond For clusters, 
the search for ultrafast phenomena presenting nonther- 
mal features is difficult. The reason for that is the lack 
of well established thermal effects from which ultrafast 
nonthermal phenomena could be distinguished. There 
is, however, a very robust thermal phenomenon in clus- 
ter physics: the cooling of vibrationally excited Cgo clus- 
ters through emission of carbon dimers 0,^. This re- 
markable process is closely related to the cage structure 
and the high stability of fuUerenes. On large time-scales, 
successive and stochastic dimer emission dominates the 
decay of a hot fuUerene. Now, a fundamental question 
regarding the short-time behavior of fuUerenes is whether 
the dimer-emission is also present in the nonthermal re- 
laxation of Ceo after an ultrafast excitation, or if other 
mechanisms become dominant. In this letter we show 
that the dimer emission is not present in the nonthermal 
response of Cgo immediately after excitation by an in- 
tense femtosecond pulse. Instead, the energy absorbed 
from the laser pulse concentrates first on a breathing 
mode and part of it is then rapidly transferred to sin- 
gle C atoms, which leave the cluster before any kind of 
thermalization of the atomic degrees of freedom can take 
place. We will refer to this process as "cold" fragment 
emission. 

Recently, several experiments on femtosecond excita- 
tion of fuUerenes have been performed [^-12|. Experi- 
mental results suggest that, as expected, the long-time 
relaxation is dominated by dimer emission. However, no 
clear conclusions can be drawn about the fragmentation 
behavior at very short times. Since such illumination by 



intense femtosecond pulses involves absorption of many 
photons by the cluster, many different electronic excita- 
tions are created, which lead to different relaxation pro- 
cesses and relaxation products. As a consequence ex- 
perimental results yield a superposition of thermal and 
nonthermal fragmentation, (multiple) ionization effects 
and even Coulomb explosion. Due to the complexity of 
the problem an exact treatment of the excitation and re- 
laxation processes is not possible. Therefore we simplify 
the problem and focus in this paper on a particular relax- 
ation mechanism, which results from excitation of many 
electron-hole pairs on the neutral cluster. This means 
that, from all possible excitations upon ultrafast absorp- 
tion of many photons we only consider those involving 
excited electrons below the continuum states. Thus, we 
do not consider single or multiple ionization processes, 
which are, of course, also present. However, from phase- 
space arguments the processes we take into account in 
this work should be the most probable ones for not very 
high laser intensities. 

In order to describe the nonthermal dynamics of Cgo we 
write down a classical Lagrangian for the atoms, which 
contains the effect of the electronic system as a many- 
body potential <I>({r.y }, i): 



N 



(1) 



Here, a number of N atoms with masses at posi- 
tions r-i interact through the potential $({r,y},t), which 
depends on the distances = jr^ — rj\ between the 
atoms. Thus, the first term of C is the kinetic en- 
ergy of the N particles, calculated from the velocities r^. 
$({ry},t) = ^{Hci{{rij}),t) is a complicated functional 
of the electronic Hamiltonian Hd- The equations of mo- 
tion for a cluster can immediately be derived through 
the Euler-Lagrange equations. In the case of the tight- 
binding (TB) Hamiltonian employed in this work the 
gradients of the interaction potential $({rij},t) are not 
available in a closed form and the equations of motion 
cannot be integrated analytically. Thus, a numerical pro- 
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cedure has to be used. We employ the Verlet algorithm 
in its velocity form [ p^p^ 

For the determination of the interaction potential 
<i>({rij},i) between the atoms we employ an Hamilto- 
nian H that consists of a tight-binding (TB) part i?TB 
for the electronic system, the other part being a repulsive 
potential (t'(j'ij) that takes care of the repulsion between 
the ionic cores: 

H=HTB+Y.<f'(^.j) (2) 

with 

Here, n^^ represents the occupation number operator for 
the orbital rj of atom i, and c^^ are fermion creation 
and annihilation operators, and the hopping matrix el- 
ement has been abridged by t^^ . For the description of 
carbon, the 2s, 2px, 2py and 2pz orbitals are taken into 
account. The angular dependence of the hopping ma- 
trix element is treated following the work of Slater and 
Koster while for the radial part of t^j and for the 
distance dependence of the repulsive potential 4>{rij ) we 
employ the form proposed by Xu et al. Diagonal- 
ization of the Hamiltonian of Eq. (||) yields the energy 
spectrum {^m{{rij{t)})} of the material and thus allows 
for the calculation of the potential energy 

= I]^(em,Oem + I]0(^y)- (4) 

Here, n{ern,t) is a time-dependent distribution of the 
electrons over the energy levels €,„. Initially, it is 
given by a Fermi-Dirac distribution nP{em) = 2/(1 -t- 
exp{(em — ^J■)/kBTc}) at a given electronic temperature 
To. Its time dependency that is caused by the absorption 
of an ultrafast laser pulse and by subsequent thermaliza- 
tion is calculated according to 

= / duj g(uj, t - At)\\n(e„i -huj.t- At) 
dt 7_oo ■ I 

+n{e,n + Tiuj,t- At) - 2n(e™, t - At)] | 

_ n{em,t) - n°(em) 
n 

(5) 

Here, the laser pulse is characterized by an intensity func- 
tion g{uj,t) which describes the distribution of intensity 
over time and energies. Thus, the electronic distribution 
TT-iem, t) is at each time step folded with the current laser 
intensity function g{uj,t). This means that at each time 
step, the occupation of an energy level Cm changes in 



proportion to the occupation difference with respect to 
levels at em — and at e^, -I- ?iu;. We model the com- 
plex processes of electron-electron collisions, that lead to 
an equilibration of the electronic system, by a rate equa- 
tion of the Boltzmann type for the distribution n(em, t). 
Thus, with a time constant ti, the distribution n{em,t) 
approaches a Fermi-Dirac distribution nP{em) at a high 
electronic temperature Tg. As we are not aware of a mea- 
sured relaxation time in Ceo, we use ti = 10 fs, a value 
that was reported for GaAs [ pT[ . 

Now we can determine the forces which are needed for 
the solution of the equations of motion by calculating the 
gradient of the time-dependent potential <I>({ry (i)}, t) of 
Eq. |: 

M{nj{t)},t) = ^Yn{ern,t){m\Vk HTBi{nj{t)})\m) 

i<j 

(6) 

where Vfc = d/dvk and \m) is the eigenvector of H corre- 
sponding to eigenvalue e™. In this equation, a term which 
follows from the gradient of the occupations n(em, t) has 
been neglected. In Eq. (^) the Hellman-Feynman theo- 
rem has been used. It is important to keep in mind that 
we are actually using a generalization of the adiabatic 
principle when we consider TB energy levels with time- 
dependent fractional occupation numbers n(em,t). Note 
that a calculation of the true nonadiabatic evolution of 
electronic wave functions is at present only possible for 
two or three degrees of freedom, while this work studies 
the time evolution of 3A^ = 180 degrees of freedom in the 
case of Ceo- 

Now we present the results for the fragmentation of Ceo 
clusters, which show important differences between ther- 
mal and nonthermal response. As experiments usually 
provide only electron-emission spectra or mass spectra 
of ionized fragments on long time scales, this theoreti- 
cal investigation is complementary to the experimental 
results in the sense that it can clarify mechanisms and 
time evolution of the damage in the clusters during the 
first stages of the relaxation process. 

A typical fragmentation process is shown in the struc- 
ture snapshots of Fig. 0, which corresponds to an ab- 
sorbed energy Eq = 3.5 eV/atom from a r = 80 fs laser 
pulse. While 40 fs after the pulse maximum, the cluster is 
still intact, already 50 fs later, the structure has been torn 
open and we can see carbon atoms and chains dangling 
from the remainders of the Ceo cage. Again 50 fs later 
the emission of three carbon monomers is observed. This 
emission of monomers is the dominant initial fragmenta- 
tion mechanism we obtain for femtosecond laser pulses. 
In the further subpicosecond dynamics of the main frag- 
ment four more monomers are emitted. Carbon atoms 
that have moved far away from the remaining cluster are 
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not shown in the subsequent panels. At t = 440 fs, a coil 
of carbon chains has formed. They stabilize to form three 
independent chains of 15, 16, and 22 atoms. This frag- 
mentation product is similar to the linked chain structure 
found as a result of thermal bond breaking The 
Ci6 cluster has a closed ring structure. Note that, apart 
from the process shown in Fig. |l|, other relaxation mech- 
anisms are possible. For instance, an absorbed energy of 
3.5 eV/atom could be enough to ionize the cluster several 
times. We stress again that we only consider excitations 
below the continuum. In the experiments, however, both 
kind of mechanisms are present. 



fs. At t = 200 fs after the laser pulse maximum the first 
breaking of bonds is taking place. This damage to the 
cluster develops over the following 400 fs into a chain of 
carbon atoms that is attached at both ends to the original 
molecule. 




FIG. 1. Nonthermal fragmentation of a Ceo in response to 
a r = 80 fs laser pulse. An energy of Eo = 3.5 eV/atom was 
absorbed. While the cluster is still intact at the time f = 40 fs 
(measured with respect to the pulse maximum), at f = 90 fs 
the cage structure has already been partially destroyed. In 
(c) we see three monomers being ejected from the cluster. A 
coil of linear carbon chains remains ai t = 440 fs after the 
emission of a total of 6 carbon monomers. 

In Fig. ^, snapshots of a Cgo cluster isomerization as 
a consequence of the absorption of Eq = 2.3 eV/atom 
from a r = 5 fs laser pulse are shown. This value of 
the absorbed energy is close to the damage threshold of 
td = 2.1 eV/atom and thus the isomerization process 
takes place on a fairly long time scale of a few hundred 




FIG. 2. Nonthermal isomerization of a Ceo in response to 
a r = 5 fs laser pulse. An energy of Eq — 2.3 eV/atom was 
absorbed. Numerous bonds are broken already at the time 
t = 210 fs after the pulse. The damage develops into a long 
carbon chain attached on both sides to the remains of the Ceo 
cage structure. 

We have analyzed the nonthermal fragmentation of 
Cgo clusters for a large range of pulse durations and ab- 
sorbed energies. We observed that a certain thermal- 
ization of the atomic degrees of freedom starts to occur 
at tth ~ 3 ps after the pulse maximum. Therefore we 
define the threshold for nonthermal fragmentation as the 
maximal absorbed energy for which no fragmentation oc- 
curs for t < tth- In Fig. |^ we show the dependence of 
the threshold for nonthermal fragmentation as a func- 
tion of the pulse duration. The fragmentation threshold 
was found to vary only slightly with pulse duration. For 
pulses of r = 5 fs to r = 300 fs duration the threshold is 
around tth = 2.1 eV/atom. We expect that after redis- 
tribution of the energy deposited in the cluster, thermal 
fragmentation may take place later even for energies be- 
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low the threshold tth- 




are approximately A = 22 A, compared to A « 7 A in the 
case of a Ceo cluster at T = 300 K. Single atomic equiva- 
lence indices with values of 330 A or more correspond to 
a dangling chain of carbon atoms that have torn them- 
selves free from the closed cage of the molecule. These 
atomic equivalence indices do not show harmonic oscilla- 
tions as most of the ai of the molecule. In the right side 
of Fig. 1^ (a), the emission of a carbon dimer C2 can be 
seen, corresponding to the two approximately parallel at 
lines with rapidly increasing magnitude. Note that the 
emission process is purely thermal; the number of elec- 
trons thermally excited above the Fermi level was always 
below 1 %. 
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FIG. 3. Nonthermal fragmentation threshold for Ceo clus- 
ters as a function of laser pulse duration. In the range from 
r = 5 fs to r = 300 fs, no clear dependence of the thresh- 
old on pulse duration is observed. A damage threshold of 
td = 2.1 ± 0.2 eV can be attributed to the entire range of 
pulse durations. 

As mentioned above, it has been found experimentally 
that carbon dimers as products of the thermal fragmenta- 
tion of Ceo dominate the fragmentation spectra Q . Cal- 
culations have been reported j|] which show that heat- 
ing of Ceo to a temperature of T = 5600 K leads to a 
fragmentation of Cgo by the emission of a dimer C2. In 
order to test the validity of the nonequilibrium fragmen- 
tation of Ceo, we thermalized Ceo clusters to a temper- 
ature T = 5600 K. We also find that the thermal frag- 
mentation process is a dimer emission on a time scale 
of several picoseconds. In order to visualize the diff'er- 
ences between thermal and nonthermal fragmentation, 
we now analyze the different trajectories with the help 
of the atomic equivalence indices These quanti- 

ties, which are instructive for the characterization of the 
vibrational excitation and structural changes in Ceo, are 
defined by 



(7) 



where ri{t) is the position of atom i. Thus, the set of co- 
ordinates of the Ceo cluster yields at each time t a set of 
60 atomic equivalence indices <Ti{t). For every atom, ai 
contains the structural information of its surroundings. 
Degeneracies of the (Ji (t) are related to the symmetry of 
the molecule. In the case of an undamaged Ceo molecule 
the high symmetry of the structure leads to a time de- 
velopment of all atomic equivalence indices in a narrow 
bundle. The atomic equivalence indices corresponding 
to a trajectory of a Ceo molecule at a temperature of 
T = 5600 K are shown in Fig. ^ (a). The amplitudes of 
the oscillations of the single atomic equivalence indices ai 
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FIG. 4. Atomic equivalence indices (Ji{t) for Cgo clusters, 
(a) Thermal fragmentation at a temperature T = 5600 K. 
The abscissa indicates the absolute time of the trajectory af- 
ter heating. At t = 6.8 ps the first irreparable damage has de- 
veloped and can be seen here as the single atomic equivalence 
indices with values ai > 330 A far above the average ai(t). 
They correspond to carbon chains dangling from the surface 
of the cluster. At t — 7920 fs the emission of a carbon dimer 
C2 can be observed as two rapidly increasing atomic equiva- 
lence indices, (b) Isomerization in response to a r = 5 fs laser 
pulse. The time evolution of the atomic equivalence indices 
corresponds to the same trajectory as the snapshots in Fig. ^ 
Immediately after the laser pulse maximum at the time t — 
the cluster expands strongly. At t = 50 fs the coherent mo- 
tion of the cluster atoms dissolves into a broad spectrum of 
individual movements. Atomic equivalence indices with val- 
ues above ai ~ 320 A correspond to chains of atoms that 
are dangling at the surface of the original cage structure, (c) 
Fragmentation in response to a r = 80 fs laser pulse. The 
ai evolution corresponds to the cluster in Fig. |l| An energy 
-Eo = 3.5 eV/atom was absorbed from the pulse. The time 
t = fs corresponds to the peak of the laser pulse. 

In order to illustrate the dramatic qualitative differ- 
ences between thermal and nonthermal fragmentation we 
show in Fig. |^ (b) atomic equivalence indices cri{t) cor- 
responding to the ultrafast isomerization of Cgo upon 
excitation with a laser pulse of r = 5 fs (see Fig. |^). 
Beginning at the peak of the very short laser pulse at 
t = fs, the cluster expands strongly. However, this 
coherent motion is quickly resolved into an incoherent 
oscillation of the cr^. Then, at lines with large values 
emerge ai t — 120 fs and they do not oscillate harmoni- 
cally. This indicates damage of the closed cage structure. 
The range of the ai values widens for subsequent times 
in accordance with the formation of a protruding carbon 
chain that was already mentioned in the description of 
Fig. |. 

Fig. H (c) shows the time development of atomic equiv- 
alence indices ai{t) during the fragmentation of Cgo (see 
Fig. During the comparatively long laser pulse of 
r = 80 fs duration, the cluster exhibits a coherent breath- 
ing mode, but shortly after the pulse maximum at i = fs 
a strong expansion of the cluster is observed. In this case 
the energy deposited in the cluster was so high that the 
disintegration sets in already ai t — 70 fs, recognizable 
by the rapidly increasing distances between the atomic 
equivalence indices ai (t) . The fact that only a small part 
of the ai lines stay relatively close to each other corre- 
sponds to the fact that a very open linear chain structure 
has formed in Fig. |l|. 

Summarizing, we have shown that the excitation of 
Cgo with femtosecond pulses gives rise to a nonthermal 
response which is qualitatively different from the well 
known thermal emission of dimcrs. Since we consider 
only electronic excitations below the continuum states, 
a comparison of our with existing experimental results 



based on detection of ionic fragments is difficult. How- 
ever, experimental fragment mass spectra [p|-p^ show 
clearly the existence of fragments of different sizes and 
a particularly large peak for C"*", which would confirm 
our simulations. In order to check that single ionization 
does not affect the main predictions of this work we have 
performed calculations on C^q clusters and we obtain es- 
sentially the same nonthermal fragmentation thresholds. 
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